Report LPSC-03-06 



Secondary proton flux induced by cosmic ray interactions with the atmosphere 
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The atmospheric secondary proton flux is studied for altitudes extending from sea level up to the 
top of atmosphere by means of a 3-dimensional Monte-Carlo simulation procedure successfully used 
previously to account for flux measurements of protons, light nuclei, and electrons-positrons below 
the geomagnetic cutoff (satellite data), and of muons and antiprotons (balloon data). The calculated 
flux are compared with the experimental measurements from sea level uo to high float ballon alti- 
tudes. The agreement between data and simulation results are very good at all altitudes, including 
the lowest ones, where the calculations become extremely sensitive to the proton production cross 
section. The results are discussed in this context. The calculations are extended to the study of 
quasi trapped particles above the atmosphere to about 5 Earth radii, for prospective purpose. 

PACS numbers: 94.30.Hn,95.85.Ry,96.40.-z,13.85.-t 



INTRODUCTION 



The study of the particle flux in the earth neighbour- 
hood has regained some interest recently with the emer- 
gence of a new generation of embarked experiments both 
at balloon and satellite altitudes. In this context, new 
measurements of the proton flux have been performed 
P, @] , providing a broad set of accurate data which could 
be compared to the latest generation of calculations. This 
is a quite compelling test which concerns the atmospheric 
flux of all secondary particles as well, since the latter are 
driven by p induced cross sections along the atmospheric 
cascade. The interaction dynamics of the incident Cos- 
mic Ray (CR) flux with the earth atmosphere and the 
earth magnetic field, for secondary particle production, 
is a complex process. So far it has been investigated 
only by means of a theoretical approach based on a dif- 
fusion equation Q . A new detailed study of this process 
through the body of recently measured data should sig- 
nificantly improve the current status of the knowledge in 
this matter, and it should validate the calculations based 
on this, or similar, simulation procedure for the evalua- 
tion of all the atmospheric secondary particle flux. It is 
also likely to improve our knowledge on the dynamics of 
the population of the radiation belts as well. 

Studying the secondary proton flux in the atmosphere 
in this context is therefore of particular interest since it 
is highly sensitive to all the components of the simula- 
tion process, in particular to the secondary proton pro- 
duction cross section as discussed below. A successful 
account of this flux through the range of atmospheric 
altitudes is then likely to cast robust grounds to this ap- 
proach in general, and to further validate the computa- 
tion techniques used. Studying this flux at higher alti- 
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tudes in the Earth environment on the basis establised 
previously should also be a useful investigation both from 
the point of view of the particle dynamics and for the fu- 
ture satellite experiments for which this background must 
be known. 

The present work is a further step of a research pro- 
gram whose previous results on the flux of secondary 
atmospheric particles at satellite altitude, in particular 
on the interpretation of the AMS 98 measurements of 
protons, leptons, and light ions, below the geomagnetic 
cutoff (GC) Jl,|10, and in the atmosphere (muons and 
neutrinos) |jj, have been reported recently. The investi- 
gation of the proton flux reported in pj is extended here 
to the atmospheric altitudes and to the high altitude re- 
gion. The article reports on the calculated atmospheric p 
flux from sea level to balloon altitude, on the comparison 
to the experimental data, and on the predictions of the 
flux at high altitudes from top of atmosphere (TOA) up 
to 3 10 4 km. The main features of the calculations are de- 
scribed in section m The production cross sections used 
in the event generator are briefly discussed in section ITTll 
The simulation results are discussed in section llVl while 
the transport equation approach is described and the re- 
sults are shown in section [3 The proton flux at high 
altitudes is reported in IVII The work is concluded in 
section TVIII 

The paper parallels a similar study on the antiproton 
flux in the earth environment 0, referred to as I in the 
following. The two papers are presented in this order for 
historical reasons. 



II. SIMULATION CONDITIONS 

The flux of secondary atmospheric protons has been 
investigated using the same simulation approach used in 
0, H la an d for the antiproton (see I) flux in the at- 
mosphere. The same computer code has been used here 
for the charged particle propagation in the terrestrial en- 
vironment including the atmosphere, as in the previous 
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studies. Incident CR proton and helium particles were 
generated and propagated inside the earth magnetic field, 
interacting with atmospheric nuclei according to their to- 
tal reaction cross section and producing secondary nuclc- 
onsp, n, light nuclei, leptons (p, e , v) from meson (n, K) 
decay, and antinucleons p , n , with cross sections and 
multiplicities as discussed below. Each secondary par- 
ticle produced in a given collision is propagated in the 
same conditions as incident CRs in the previous step, 
resulting in a more or less extended reaction cascade de- 
veloping through the atmosphere, which included up to 
about ten generations of secondaries for the protons of 
the simulation sample Q. 

The reaction products are counted each time they are 
crossing, upwards or downwards, the virtual detection 
spheres. The locations of the latter were chosen between 
sea level and about 36 km for ground and balloon ex- 
periments (BESS, CAPRICE), at 380 km for the AMS 
satellite experiment, and beyond up to about 3 10 4 km for 
the rest of the study. All charged particles undergo en- 
ergy loss by ionisation in the atmospheric medium. Each 
event is propagated until the particle disappears by nu- 
clear collision, stopping in the atmosphere by energy loss, 
or escaping to the outer space beyond twice the genera- 
tion altitude (see [4j,|5j,|6J and I). 

The incident CR proton and helium flux have been 
measured recently by several experiments 0, 0, 0, 0, 
H^ |. In the present work, functional form fits to the AMS 
data 0, |R| have been used in the calculations to gen- 
erate the corresponding flux. For other periods of the 
solar cycle than those of the measurements, the incident 
cosmic flux were corrected for the different solar mod- 
ulation effects using a simple force law approximation. 
The A > 4 components of the CR flux were not taken 
into account in the present calculations (see ref Q for 
details). 



IV. RESULTS 
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FIG. 1: Simulated proton flux in the atmosphere (histograms) 
compared to experimental results (symbols) between TO A and 
sea level. From top to bottom: CAPRICE TOA @/ (fM 
squares), CAPRICE 94 Q/ (7 altitudes: 29. £22.1, 17.5, 16, 
12.6, 9.9, 5.75, in km), Kocharian et al, Jwj (3.2 km, full 
triangles), BESS [2] (2.77km, full squares), Diqqory et al, 
Jfl'j] (sea level, inverted full triangles). See also \ll\l for other 
data at sea level. 



III. CROSS SECTIONS 

The inclusive p + A — > p + X and 4 He + A — > p + X 
proton production cross sections and the proton and 4 He 
total reaction cross sections on nuclei used here are de- 
scribed in I. For incident protons however, the description 
of the production cross section below 7.5 GeV [14j used in 
has been constrained by the 4.2 GeV measurements 
from ref [l^. This decreased the cross section over this 
range of energy by a few tens of percents. The details 
will be reported later. 

The neutron production cross sections were taken the 
same as for protons. This is expected to be a good ap- 
proximation for the incident and final state energies con- 
sidered here. Charge-exchange reaction channels have 
been neglected in account of their small cross sections. 
The nucleon production from He fragmentation were not 
included. 



The proton flux in the atmosphere have been mea- 
sured recently by the CAPRICE experiment between 5 
and 29.9 km Q and at TOA Q, and by BESS at lower 
altitude (2.77 km) 0, while previous measurements were 
available from 0] at 3.2 km, and from a t sea level. 
This collection of data points is compared on figurenwith 
the simulation results between sea level up to balloon 
altitudes. The comparison is remarkably good through 
the whole range of altitudes. Note that the calculations 
become increasingly sensitive to the proton production 
cross-section when going from high to low altitudes since 
all secondary protons result on the average from a se- 
quence of n collisions N + A — > p + X, N standing for 
nucleon and n being the average collision rank, involving 
(approximately) the same inclusive nucleon production 
cross section for each collision. Let the latter be noted 
ag dcr N (Eu8,E) a gj ven co nj s j on w jth Ei, E, 9, be- 
ing the nucleon incident energy, proton final energy, and 
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FIG. 2: Rank distributions of the detected protons obtained from the simulation, for a set of altitudes between sea level ana 
balloon high float altitude, showing the evolution of the mean production rank in the collision sequence, of the detected protons. 



production angle, respectively. Although the collisions 
involve different values of the kinematic variables, the fi- 
nal proton flux for a given altitude depends qualitatively 
on some appropriate average of the cross section above 
to the nth power, i.e.,: < d<TN ^2^' E ^ >™j with n increas- 
ing with the decreasing altitude. Therefore the lower the 
altitude, the more sensitive the flux is expected to be to 
the inclusive proton production cross section. 

Fig [21 shows a sample of rank distributions from sea 
level to ballon altitude obtaines from the simulation sam- 
ple. The mean rank n appears to drop from around 7 at 
sea level down to about 1.5 at 36.5 km. The flux is thus 
sensitive to approximately the 7th power of the cross sec- 
tion at sea level, while it depends only to the 3/2 power 
of this cross section, approximately, at ballon altitude. A 
discrepancy of a factor of 2 between calculation results 
and data at sea level can thus be considered as a good 
result since it would point to a discrepancy on the cross 
section of less than 10 percent. 

These calculations thus provide an extremely sensitive 
test of the correctness of the cross sections used, and be- 
yond of the overall method, and they assign a high con- 
fidence level to the all secondary particle flux calculated 



by this approach. 

V. THE TRANSPORT EQUATION APPROACH 

The proton flux has also been calculated u sing the 
same diffusion equation approach as used in refs [n| (see 
also ) in order to compare the two methods. 

A. Numerical integration method 

The equation to be solved can be written as: 

dNi,_d_/dEA Mt jp \ , N t (x,E t ) 
dx + dEi\dx Xmt {E t ) 
-EaJ^-J $E th m-(E l ,E A )N A (x, E A )dE A = (1) 

Where i represents the transported particle species, 
Ni(x,Ei) its energy (Ei) dependent flux after crossing 
the thickness between x and x + dx (in g/cm 2 )of (atmo- 
spheric) matter. The second term is the particle energy 
changing term accounting for energy loss by ionization 
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(using the Bethe-Bloch formula). In the third (absorp- 
tion) term, X int is the particle interaction length derived 
from the total reaction cross section of the considered 
system. The last term is the source term, accounting 
for particle creation, with A standing for the projectile- 
target system leading to the production of particle i. 
4§:(Ei, Ea) is the production cross section for particle 
i in the system A with incident energy Ea- The sum 
runs over all the allowed A channels. The differential 
cross sections used here are the same (angle integrated) 
as used in the simulation. (m a i r ) is the mean nuclear 
mass of the atmospheric nuclei (14.58 amu). 

The equation has been solved for the secondary proton 
flux at various altitudes. The technique used proceeds 
by the method of finite differences with implicit scheme 
proposed in ref. . It consists of solving a system of 
linear equations: 
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Where m and n are the index of the steps in energy and 
thickness crossed respectively. 

For practical reasons related to the energy range to 
be covered versus the number of steps required and the 
size of the matrix to be inverted, the definition of the 
derivative used at the boundaries had to be modified to 
obtain stable calculation results 221. 
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FIG. 3: Results obtained using the diffusion equation approach 
for the transport of the proton flux in the atmosphere. The 
data are the same as in Fig. Q 



B. Results 



Fig.[3]shows the results obtained in these calculations, 
compared to the same experimental data as in Fig. ^ 
at various altitudes in the atmosphere. The agreement is 
seen to be fair in the high altitude range, down to around 
10 km. For the 5.5 km data, the observed disagreement 
is similar to that obtained by simulation, the data points 
having a much different energy behaviour than observed 
for the higher and lower energy spectra. 

For altitudes lower than 5 km, significant disagree- 
ments between data and calculations appear, the latter 
overestimating the data by a factor of about 2 to 4. This 
disagrement probably originates in the one-dimension ap- 
proximation of this approach. The difference between the 
3-dimensional (3-D) approach of the simulation and the 
1-dimcnsion diffusion equation lies in a few 3-D effects 
not included in the latter approximation : a) The curva- 
ture of the particle in the magnetic field, b) The angle de- 
pendence of the particle production, 1-D approximation 
using the angle integrated cross section. A consequence 
of the above is that the effective path of particles in the 



atmosphere is longer than assumed in the 1-D approxi- 
mation. Using the more realistic values of the path ob- 
tained in the simulation, in the diffusion equation brings 
only minor improvements to the observed disagreement. 
Introducing some cuts in the angular range of integration 
of the cross section distorts the resulting spectrum and 
does not improve the discrepancy either. 

It is interesting to note that the failure of the 1-D ap- 
proximation reported above at low atmospheric altitudes 
seems not to have a significant effect on the goodness of 
this approximation for the calculation of the atmospheric 
neutrino flux. This conclusion was reached on the basis 
of the fair agreement obtained between the 3-D and 1-D 
calculations in Ref.0. This can be understood by the 
fact that the low altitude neutrino production represents 
less than 10% of the whole atmospheric production (frac- 
tion decreasing with the decreasing neutrino energy). 
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FIG. 4: Simulation results for the altitude dependence of the downgoing (left) and upgoing (right) proton flux. Each set of 
curves of the two panels corresponds to a bin of geomagnetic latitude 8 m as indicated on the figure. The points and curves show 
the energy integrated primary flux of protons (full triangles), and the flux for proton with kinetic energies above 0.1 GeV (solid 
line), 0.3 GeV (dashed line), 1 GeV (solid line), and 3 GeV (dash-dotted line), respectively. See text for the discussion. 



VI. SECONDARY PROTON FLUX VERSUS 
ALTITUDE 

The same study of the particle flux at high altitudes as 
performed in I for antiprotons has been conducted here 
for protons over a range of altitudes going from TOA up 
to around 3 10 4 km, with the same purpose of, on one 
hand, understanding better the dynamics of the proton 
population at high altitudes, and on the other hand to 
provide reliable predicions of the atmospheric secondary 
proton flux for future embarked experiments. 

Note that the calculations reported in the following 
include only secondary protons produced in the atmo- 
sphere. They do not include the contribution originating 
from the atmospheric production of neutrons decaying 
into protons. This contribution is negligible inside the 
atmosphere, but it is a well known major component of 
the radiation belts for the altitudes considered below. 
This latter issue is investigated in a separate study. 

Fig. 2| left shows the angle (8 Z < ir/4 rad) and energy 
(E > 0.1 GeV) integrated downgoing primary proton flux 



from sea level up to 3.2 10 km, in bins of geomagnetic 
latitudes between equator and polar region. At low lati- 
tudes, the energy integrated flux is observed to drop by 
more than one order of magnitude between asymptotic 
distance from Earth and about 10 3 km because of the 
geomagnetic cutoff which is most effective at these lat- 
itudes. In the polar region where there is no GC, the 
flux is predicted constant over the same range of altitude 
until it reachs the atmosphere. Inside the atmosphere 
it drops exponentially due to particle absorption. The 
secondary flux distributions are shown on the same fig- 
ure for different energy integration thresholds. The main 
features of these distributions clearly do not depend crit- 
ically on the energy threshold (0.1 to 3 GeV), although 
significant differences are observed, discussed in the fol- 
lowing. The broad peak at low altitude corresponds to 
the atmospheric production yield, the low altitude side of 
the peak being governed by atmospheric absorption, and 
the high altitude side by the atmospheric density. The 
intermediate plateau of the flux distribution corresponds 
to the population of quasi-trapped particles which are 
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FIG. 5: Same as figure |7| for the energy integrated downward 
proton flux (full triangle) and the E> 0.1 GeV (full circles and 
solid line). The full circles symbols and dashed line show the 
fraction of the latter for which the first adiabatic invariant 
conservation is violated. 



accomplishing a few rides between mirror points before 
being absorbed in atmosphere (see discussion in I). The 
peak observed at high altitude around 20-30 10 3 km for 
the low latitude region, corresponds to low energy par- 
ticles for which the first adiabatic invariant is violated 
(see Fig. |5}, and thus drifting to higher shells without 
being rapidly absorbed in the atmosphere on the normal 
trajectory between mirror points (see examples in I). All 
distributions have a high altitude cutoff, which is both 
energy and latitude dependent, dropping from 3 10 4 km 
for low energy low latitude particles, down to 10 2 -10 3 km 
at high energy and high latitude. This latter feature can 
be understood qualitatively by simply considering that: 

a) Particles produced at low latitudes have a natural mo- 
mentum limitation set by the simple condition that their 
gyration radius is smaller than the distance of the gyra- 
tion center (mean field line) to the atmosphere. 

b) Particles produced in the polar region will tend to es- 
cape in account of the low value of the magnetic field at 
the poles, and subsequently, the higher the energy, the 



more effective the trend. 

Quantitatively, this high energy cutoff can be understood 
in the Stormer approach to the problem |23| which pro- 
vides the maximum energy allowed for a trapped particle 
at given altitude and latitude. For trapped particle, the 
predicted dependence of this GC on the distance, on the 
energy, and on the latitude, completely accounts for the 
features observed in Fig. 21 

Integrating the zenith angles over a larger range of ac- 
ceptance produces distributions significantly wider, ex- 
tending to lower altitudes, as it could be expected from 
the above considerations. 

This is further illustrated on Fig. fright which shows 
the same distributions for the upgoing proton flux, com- 
plementing the previous figure. The upward primary flux 
begins around 300 km in altitude at the upper bound- 
ary of the forbidden region where primary trajectories 
are not allowed by the GC 23]. It converges asymp- 
totically with the incoming (downward) primary flux at 
the highest altitude calculated (primary flux isotropy). 
For the upgoing secondaries, a similar plateau and high 
altitude peak at low latitudes, are observed as for the 
incoming secondary flux. These parts of the flux distri- 
butions closely overlap with those of the incoming flux 
at low and intermediate latitudes and for the lower en- 
ergy particles, showing that it is a flux of quasi-trapped 
particles. At high latitude the high energy flux which is 
limited to the inner atmosphere for downgoing particles, 
is to the opposite located at high altitude above 100 km, 
i.e., approximately above TOA, for protons above 3 GeV. 
These are escape particles which can be produced only 
tangentially to the atmosphere - since forward produced, 
the backward production cross section being vanishingly 
small - to have a chance of being deflected upwards by 
the magnetic field (with a strong east-west asymmetry 
effect, see p|). Consequently, the particle trajectories of 
this flux have a large zenith angle on the average. 

Fig. [5] shows three distributions with the same lati- 
tude binning as on the previous figures: The angle and 
energy integrated downgoing proton flux from sea level 
up to 3.2 10 4 km (full triangles), the energy integrated 
(E > 0.1 GeV) secondary proton flux (full circles, solid 
line), these two distributions already seen before, and 
the fraction of the secondary flux for which the conser- 
vation of the first adiabatic invariant (magnetic momen- 
tum of the particle (2^) is violated by more than one 
order of magnitude. For the low latitudes and high alti- 
tudes above about 10 3 km, the secondary flux is almost 
exclusively adiabatic invariant violating, while the con- 
servation of the invariant is approximately satisfied for 
quasi-trapped particles over the range of altitudes be- 
tween I0 2 km and I0 3 km. 



VII. SUMMARY AND CONCLUSION 

In summary, it has been shown that the simulation ap- 
proach to the proton flux in the atmosphere allows to sue- 
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cessfully reproduce the data to a high level of accuracy. 
This result confirms the reliability of the method. The 
proton flux has been calculated up to around 10 Earth 
radii. The results have shown that a large component of 
quasi-trapped particles dominates the flux over the inter- 



mediate range of altitudes (10 2 -10 4 km). These results 
should serve as a guideline for the evaluation of the par- 
ticle background for experiments embarked on satellites. 

A study of trapped protons originating from the neu- 
tron flux is in progress, in the same context. 
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